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We consider the spin physics of fractional quantum Hall effect as a continuous function of the filling 
factor, and predict a rich phase diagram with magnetic phase transitions whose phase boundaries cut 
across many fractional plateaus. The calculated phase diagram also applies to fractional quantum 
Hall effect in systems where valleys may play the role of spin, such as AlAs quantum wells, graphene, 
and H terminated Si(lll) surface. 
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The interplay between the Coulomb interaction and 
the electron spin degree of freedom has led to remark- 
ably rich physics in the fractional quantum Hall efFeclfS 
(FQHE). Phase transitions have been observed as a 
function of the Zeeman energy, E-z, in transport,'2Hsl 
opticaPltiSl and NMR experiments,'^^!^^^ establishing that, 
in general, FQHE states with several discrete values of 
polarization can occur at a given filling factor. The ba- 
sic physics is understood. The FQHE at filling factor 
V = n/ (2rt± 1) is a manifestation of the integer quantum 
Hall effect (IQHE) of composite fermions at filling n. For 
n > 1, several states with n — -\-n^ are available, with 
polarization P = {n^ — n^)/{n^ + n^). Which of these 
states is favored as a function of Ei has been determined 
in detailed microscopic calculations within the composite 
fermion (CF) theory,'iSHlfila,nd can qualitatively be under- 
stood in terms of a competition between the Zeeman and 
the CF cyclotron energies. 

The above physics carries over to FQHE in systems 
with two valleys (with the actual spin frozen due to large 
Zeeman energy) , with the valleys playing the role of spin 
and the valley splitting of i?z- Transitions between differ- 
ently valley polarized CF states have been seen in AlAs 
quantum wells,!^ where the valley splitting is controlled 
by application of an in-plane strain. In high quality H- 
terminatcd Si(lll) surface with two degenerate valleys, 
the 2 — rt/ {2n ± 1) states with n ~ 2, 4j6 have been ob- 
served but those with odd n are absent.'^ This behavior 
is consistent with what is expected in the limit of vanish- 
ing valley splittin^^^ if one neglects the residual interac- 
tion between composite fermions, a gapped state occurs 
only for even n. In graphene as well, an effectively two- 
component CF physics is observed in certain regimes; e.g. 
the experiments of Refs . 1^ and show that in the range 
1 < < 2 FQHE occurs at 2 - n/{2n ± 1) for n = 2,4 
but not for odd n (although both even and odd n are 
observed for i/ < 1). Our calculations and results below 
apply to FQHE in any system where two components 
are active, but we will use, for simplicity, the convenient 
language of "spin" and "Zeeman energy." 

Theoretically, to date, the SU(2) physics of only the 
FQHE states at the discrete filling factors n/ (2n± 1) has 
been studied.'^^HUl jj^ this Letter we obtain the phase di- 
agram as a continuous function of the filling factor. This 
brings out the fuller richness of the spin physics and re- 
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FIG. 1. Schematic description of the states with different 
polarizations in the filling factor range 1 > i^* > 2 (top), 
2 > V* > 3 (middle) and 3 > v* > 4 (bottom). Composite 
fermions are pictured as electrons carrying two arrows (rep- 
resenting bound vortices). The horizontal lines depict the CF 
ALs, with the spin-up ALs on the left and spin-down ALs on 
the right. The states are labeled (z/^,!/|); the filing factor of 
the partially filled A level is 



veals additional phase transitions. Further investigation 
of this phase diagram and the nature of the associated 
phase transitions should produce deeper insight into the 
underlying physics of the FQHE. 

At an arbitrary filling in the range 1/3 < v < 2/3, 
electrons ca pture two vortices to transform into compos- 
ite fermions j'^SEH which form Landau-like levels called 
A levels (ALs), and have a filling factor i/* given by 
1/ = V* l(2v* ± 1). We allow occupation of spin-up 
and spin-down ALs, with v* = ~\- v*^. The corre- 
sponding state, labeled (y^,v*^)^ has spin polarization 
P = (t/| — v^)/!^*. (We take below the convention in 
which the up spin is favored by the Zeeman energy, i.e., 

^ ^l') We will make the assumption, valid for weakly 
interacting composite fermions, that only one of the ALs 
is partially occupied, so one of vi- and i/f is an integer 
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FIG. 2. The total energies of various states in the 

filhng factor range Xj'i < v < 4/9. The energy E is quoted 
in units of /el, where I = ^Jhc/eB is the magnetic length 
and e is the dielectric constant of the background material. 
The dark (open) symbols correspond to states in which the 
partially filled AL contains a crystal of CF particles (holes). 
In this and subsequent figures, the scale for filling factor is 
different in the three columnar regions. 

and the other will be written as j ' + P, where v is the 
filling factor of the partially filled AL. Fig. [T] depicts the 
possible states between 1/3 and 4/9 studied below. 
The wave function of the state is given bj^^ 

= ^.{zi, • • • zn}ui ■ ■ ■ UN^dN^+i ■■■(In] (1) 

N 

-^LLL n (2) 

j<k=l 

= '^^'^{zir ■ ■ zn^}^ui{zn^+i - ■ ■ zn} (3) 

Here Zj — Xj + iyj denotes the position of an electron; 

and $1^1 are wave functions of spin up and spin down 
fermions at and v'^] u and d are the up and down 
spin spinors; A indicates antisymmetrization; and Plll 
denotes lowest Landau level (LLL) projection. For the 
states relevant to our current study (Fig. [T]), all occupied 
orbitals of spin-down electrons in ^^■^ are definitely occu- 
pied for spin-up electrons ; the product wave function 
is therefore annihilated by the spin lowering oper- 
ator (i.e., it satisfies the Fock conditiorP^l) and is thus an 
eigenstate of both Sz and S"^ with S = Sz- These proper- 
ties are preservecP^ under composite fermionization via 
Eqs. 1 and 2. The LLL projection will be carried out 
using the method described elsewhere in the literature.'^ 
The interaction energy 

where V is the interaction, will be evaluated by the 
method of Metropolis Monte Carlo. 



FIG. 3. Difference between the energies of the states with dif- 
ferent polarizations, shown near each line, in the filling factor 
range 1/3 < u < 4/9. 



Before proceeding further it is necessary to know the 
state of composite fermions in the partially filled AL (as 
is the case for one of the two spins), the nature of which 
is dictated by the weak residual interaction between com- 
posite fermions. Close to integer values of v* , we expect 
the formation of a crystal of CF particles or holes. We 
will assume that the CF crystal is an adequate model for 
the entire range < P < 1, and for completeness, we 
will consider both the CF particle and CF hole crystals 
for the entire range. It is possible that some other state 
of composite fermions (a FQHE liquid, a stripe phase, 
a paired sate, or a crystar^) may have a slightly lower 
energy than the CF crystal at some intermediate fillings, 
but we will argue below that the energy differences be- 
tween such states at a given spin polarization are small 
compared to the energy differences between differently 
polarized states, and therefore can be neglected when our 
primary focus is on the spin polarization of the state. 

We will use the spherical geometrjl^H for our calcula- 
tion, in which N electrons are confined to the surface of 
a sphere with 2Q flux quanta passing through it. This 
geometry is convenient as it has no complications due 
to edges and enables a simple treatment of the inter- 
action with the background charge (taken as a single 
positive charge at the center). The crystal of compos- 
ite fermions is obtained by placing them at the positions 
that minimize the Coulomb energy of classical charged 
particles, according to J.J. Thomson's famous plum pud- 
ding model of the atom (these positions are tabulated in 
the literatur^2ll) ; as one may expect, this is essentially a 
triangular crystal with a few defects introduced by the 
curvature of the spherical surface. 

The technical details are as follows. The wave func- 
tions for the CF particle or hole crystal a.t i/* = j + D is 
constructed from the analogous electron or hole crystal 
a.t ly = j + D. We indicate the electron coordinates var- 
iously as r = {0,(j)) = {u,v), where the spinor variables 
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are defined as u = cos(6•/2)e''^/^ i; = sm{e /2)e-"t'^^ . 
For the electron crystal, we construct the Slater determi- 
nant wave function with j Landau levels (LLs) 
completely occupied, and electrons in the partially filled 
+ LL localized at crystal locations. The wave func- 
tion of an electron wave packet in the (j -I- 1)*^ LL (where 
j = is the LLL) localized at {U, V) is given by^SI 

/,2Q - (Q-\jj,^(Q+i) 



V(u.y),j - ' V{u.y),o 



(5) 



where i/'^^y) o = (^*" + is the wave func- 

tion of an electron wave packet localized at ([/, V) 
in the LLL, and S~ = v*d/du — u*d/dv is the 
LL raising operator (which also lowers the monopole 
strength by a single unit). For the hole crystal, 
we first construct a Slater determinant wave function 
<i>j+i(ri, r2, • • ■ t^m-ItTm) of j -f- 1 filled LLs. Creation 

of M — TV holes localized at • • • ,Rm-n is accom- 
plished by the action of ^(i?i) • • • ^{Rm-n), which pro- 



duces $j+i(ri,r2, ■ ■ ■ ,rN,Ri,- 



), i.e. M - N 



electron coordinates are replaced by variables denoting 
the hole positions. This wave function, however, also pro- 
duces holes in the lower LLs. That can be taken care of 
by replacing in the Slater determinant the entries x(-R) 
corresponding to lower LLs by zero, which guarantees 
that the wave function has holes only in the topmost 
LL. For example, the wave function $2 with two holes at 
= {Ui,Vi) and R2 = {U2,V2) in the 2"^^ LL is given 

by: 
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(6) 

and sub- 



With these wave functions we construct ^u'^i^i , 

While the full form of the 
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stitute in Eq. 2 to get v]/ 

LLL projected ^E*^^ is rather complex, it is well defined 
and produces accurate results. 

Fig. [2] shows the Coulomb energies of the various 
states as a function of the filling factor in the range 
1/3 < V < 4/9. The filled (open) symbols represent the 
energies for states where the partially filled AL is mod- 
eled as a crystal of CF particles (holes). (In the range 
3/7 > V > 4/9 only the CF particle crystal has been eval- 
uated.) The energies are shown for a system with = 84 
particles, but we have confirmed, by studying various A^ 
values, that the A^ = 84 results are close to the ther- 
modynamic limit. The energies include the interaction 
with the uniform positive background. The energy dif- 
ference between the CF particle and hole crystals for a 
given polarization P is small compared to that between 
states with different polarizations. Fig. [3] depicts the 
energy differences between two consecutive states, where 
we take the lower of the particle or hole crystal energies 
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FIG. 4. The calculated polarization phase diagram of vari- 
ous states iy^jU'^) in the CF filling factor range 1 < u* < A 
(1/3 < < 4/9) as a function of «: = E-Ljif? jet). The thick 
solid fines are for a system of zero thickness; the thin dashed 
lines correspond to a quantum wefi width of 40 nm for density 
lO^^cm"^; and thin green lines indicate the phase boundaries 
predicated by a model of noninteracting composite fermions. 
The shaded regions schematicafiy depict filfings where tran- 
sitions take place from one plateau to the next. The inset 
shows the polarization P = (z/^ — I'D/v* as a function of v* 
for fixed values of k = 0.0225 (black solid line), 0.0125 (purple 
crosses) and 0.005 (red circles), as evaluated for a zero thick- 
ness sample without disorder. The blue triangles indicate the 
spin phase transitions at i^* = n identified previously.'^ 



for a given spin polarization. From the difference, it is 
straightforward to determine the critical Ez = K{e^/e£) 
for the various transitions, plotted in Fig. |4j (In both 
Figs. [3] and |4] the curves have been smoothened to elim- 
inate fluctuations due to finite system size. The uncer- 
tainty in the curves due to finite size effects is estimated 
to be Sk « 0.001.) The phase diagram in Fig. |4]is the 
main result of our paper. 

The most striking feature of the phase diagram is the 
presence of lines that cut across many FQHE plateaus. 
A free-CF model provides useful insight into the struc- 
ture of the phase diagram. This model considers non- 
interacting composite fermions, with the CF cyclotron 

at 



energy taken aPI^Sl Sw! 



m* 2i/*±l — a{2u*±l) el 

where uJc — eB/nieC is the cyclotron en- 



2i/'±l' 

ergy, TOg is the electron mass in vacuum, and m* is the 
CF "polarization mass.'^^E^IIZI (Por CaAs parameters, we 
have TO*/me = 0.026a-y/i3|Tj.) This model predicts that 
the phase boundaries between the different polarizations 
are simply given by Ez = integer x huj* : 



integer 
a{2iy* ± 1) 



(free CF model) 



(7) 



where the integer goes from 1 to mt{i^*). We show the 
phase diagram predicted by the free-CF model in Fig. |4] 
(thin lines), fixing a = 15.9 by requiring coincidence at 
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the lower transition at = 2/5. (The polarization mass is 
known to be much larger than the "a.c tivation mass" rele- 
vant for transport experiments! '' I ^^ l ^ The free CF model 
reproduces the qualitative features of the phase diagram 
obtained from the full microscopic theory; the quantita- 
tive deviations between the two can be attributed to the 
residual interaction between composite fermions. 

We have also studied the effect of finite thickness as- 
suming a cosine shaped wave function in a square quan- 
tum well, which leads to a width dependent effective in- 
teraction. As an illustration, the phase diagram is also 
shown (dashed lines) for a quantum well of width 40 nm 
at a density of 10^^ cm~^. Lower densities and shorter 
widths produce a smaller deviation from the zero width 
results. 

Disorder, not included above, will affect our results in 
several ways. Most importantly, the first order phase 
transitions at the phase boundaries will turn into con- 
tinuous percolation transitions in the presence of disor- 
der, producing extended states that allow identification 
of the phase transition in transport experiments (below). 
Disorder will also affect the energies of the various states 
differently, and thereby modify the phase boundaries. Fi- 
nally, disorder will create spatial variations in the filling 
factor, which will provide a correction to the "ideal" po- 
larization. For example, while the ideal 1/3 state is fully 
spin polarized for all k, disorder will slightly diminish P 
for K < 0.03. A reliable quantitative estimation of these 
corrections is not possible at this stage. 

For sufficiently large k, the state is fully polarized at 
all fillings. Interestingly, for a range of k values, the 
state at two consecutive fractions along the sequences 
n/{2n + 1) is fully spin polarized but partially polarized 
at intermediate fillings. The measurements of Tiemann 
et al. clearly show such behavior where both 1/3 and 
2/5 are fully spin polarized but the intermediate state is 
partially polarized. The inset in Fig. |4] shows the spin 
polarization P = {ly^ — i>^)/{h'^ + '^l) ^ ^ function of v 
for several fixed values k (as would be the case when v is 
varied by changing the density at a fixed B), indicating 
complex K dependent behavior. 

Transport measurements can identify the phase bound- 
ary through the appearance of an R^x peak as a function 
of the Ez-, which appears because of the presence of ex- 
tended states (due to disorder) at the transition point. 
These R^x peaks are flanked by two states with the same 



quantized Hall resistance Rxy- (These should be dis- 
tinguished from the peaks atP^ i/* = n -|- 1/2, shown 
schematically as vertical shaded regions in Fig. |4j which 
indicate transition between states with different values 
of the Hall resistance.) Such peaks have been seen in 
a number of experiments, both at the s pecial fr actions 
ly = n/{2n ±1) and slightly away from itPSEEH Krauss 
et ali^ have detected a sharp Rxx peak inside the resis- 
tance minima near 2/3 and 3/5, and determined a phase 
boundary of the type calculated here, clearly delineat- 
ing states with various spin polarizations as a function of 
the filling factor. Unfortunately, a technical difficulty in 
dealing with reverse flux attachment makes the calcula- 
tion of the phase diagr am al ong the sequence n/{2n — 1) 
much more challenging j^^l^^ l and outside the scope of the 
present work. However, the phase transition at 2/3 is 
found to occur in Ref . [T3l at ~ lOT, which, according to 
the free CF model, corresponds to a ^ 19, which is in 
the same ball park as the value obtained above (and, in 
fact, somewhat higher, as expected from finite width cor- 
rections). Resistance spikes have also been observed at 
magnetic phase transitions in the integer quantum Hall 
regime.!^ We also note that because the phase bound- 
ary involves a change of magnetization, it should ex hibit 
hysteretic behavior, as observed pr e viously! ^ * ^ * ^ ^ 

Before ending, we note that we have focused here on 
magnetic transitions, and not attempted to discriminate 
between various possible states of composite fermions 
with the same spin polarization, such as their FQHE, 
paired state, stripes or crystals. The difference between 
the energies of such states for a given spin polarization is 
likely to be much less than that between the states of dif- 
ferent spin polarizations, as evidenced by the relatively 
small difference between the energies of the CF parti- 
cle and hole crystals. A more detailed phase diagram 
will contain additional phase boundaries separating dif- 
ferent states with the same spin polarization. Further, 
for graphene and other multi-valley systems, it will be 
important to consider bo th the spin and valley indices to 
bring out the full physics.^^ilSll 
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